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2.7 Aerosols and coagulation 



[Refs]: 



Present, Kinetic Theory of Gases 
Puchs, Mechanics of Aerosols 

Friedlander, Smoke, Dust and Haze 

Seinfeld, Atmospheric Chemistry and Physics of Air Pollution 
Levich: Physio-Chemical Hydrodynamics 

2.7.1 Brownian diffusion of particles 

Particles of sub-micrometer size collide with air molecules randomly, and behave collectively 
as a gas. 

Let n{z) — number of particles per unit volume, i.e, the number density of Brownian 
particles in air. By the perfect gas law. 



L — Avogadro's number — 6.025 x lO^^molecules /mole = 2calories /degree /mole 
When the cloud of particles is in hydrostatic equilibrium. 



p = nkT, 



(2.7.1) 



where k is Boltzmann's constant 



k = R/L, 

R — Universal gas constant = 8.317 x lO^ergs/degree/mole 



and 



dp 

dz 



-pg = -nmg 



(2.7.2) 



where m is the mass per particle. Combining the preceding two equations. 



dp mg 



p kT 



hence. 




(2.7.3) 



and 




(2.7.4) 
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which is Bolzmann's law. 

Alternatively, random collisions on the microscale give rise to diffusion on the macroscale. 
At the equilibrium state, diffusion and gravitational convection must balance each other so 
that, 

-D^ -nV^O (2.7.5) 
dz 

where V — fall velocity. Equating Stokes drag with the particle weight QTr/iaV — mg, we get 



V = ^ (2.7.6) 



Solving (2.7.5) and using (2.7.6), 



Upon comparison with (2.7.4), the Brownian diffusivity can be identified 

kT 

D^f^ (2.7.8) 

This formula, due to Einstein (1905) and Smoluchowski (1906), is valid if 2a is smaller 
than the mean free path i of air molecules. Otherwise Cunningham's empirical correction is 
needed. 



where 



a 



( \.\a\ 
1.257 + 0.4 exp f 



(2.7.10) 



is the correction factor. 

For aerosol particles in air under normal temperature, the diffusivity is 



^ 1 X 10-3 10-2 0.325 X 10-^^ , „ , 

(2.7.11) 

2a 0.15 X 10-3 a ^ ' 

As a rough order-estimate (Levich), we use Eistein's formula for water at room temperature, 
the Brownian diffusivity for colloidal particles is 

^ 0.55x10-13 

(2.7.12) 

a 

In the following table taken from Seinfeld, p. 325, D and v are compared. For gases, the 
mean free path is typically £ = 10~^ ~ 10-^ cm. 
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2 a (jU. m) 


D (cmVs (T = 20°) 


V (cm/sec) (p^ = Ig/crn?) 


0.001 
0.01 
0.1 
1 

10 

100 


5.14 xlO-2 
5.25 xlO~^ 
6.75 xlO-6 
2.77 xlO-^ 


8.62 xlO-5 
3.52 xlO-3 
3.07 xlO-=^ 
30.3 



2.7.2 Coagulation due to Brownian diffusion 

When small particles are bounced around randomly by surrounding fluid molecules, they 
may come so close to one another that Van der Waals force binds them together. This is 
coagulation. In a moving fluid additional factors such as fluid shear and Columb forces may 
intervene. A simple model (by Smoluchowski) for a stationary fluid with identical spherical 
particles of radius a is as follows . 




Figure 2.7.1: Left: A spherical shell. Right: Two spherical particles in collision. 



Let us focus attention on a flxed particle. Consider a spherical shell from r to r + rfr. 
Figure 2.7.2-left. The rate of increase of particles inside the shell is 

dt 

which must be equal to the net influx through the two surfaces of the shell 

-— -A-nr^D— dr 
or \ or 



thus. 



dn D d f 2^'^\ f2 7 13) 

dt dr \ dr I 
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Assume that whenever two particles come into contact they stick to each other and become 
one. Therefore the spherical surface of radius 2a concentric with the stationary particle acts 
as a sink, on which n — 0, i.e., 

n = r = 2a 
n — > rioo r — > oo 
See Figure 2. 7. 2- right. The initial condition is 

^ = '^oo-i 2a < r < oo t = 
Th solution can be facihtated by introducing 



Uoo-n / r 



w 



2a 



and 

it is shown in Appendix A that 
where 

The boundary conditions become, 
while 

The initial condition is 



r-2a 



X — 



dw 
'di 



2a 



D' 



D 



{2ay 



w — 1, X = 0, 

w = 0, X = oo 



u; = 0, t = 0,x>0 
The solution, which can be obtained by the similarity method (see Appendix B) 



w = 1 — erf 



X 



,2V Dt 



2 / 2VDi _z2 

e dz 



12 ji 

TT Jo 



or. 



or 



rioo V TT yo 



n 



noo 



2a 



'2 

TT Jo 
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Finally, the number concentration near a fixed particle is 



n{r, t) 



= 1 



2a 2a 
— + — 



'2 

TT Jo 



2a 2a ^ 
1 erf 



r — 2a 
2Dt 



(2.7.25) 



We now use this information to find the rate of coagulation when all particles are moving, 
by calculating the number density of particles in the process of collision. Starting from one 
particle, the rate of flux of particle across the sphere of radius r = 2a is 



J{t) = Ai^r'^D 



dn 
dr 



When 



r=2a 



t > 



AT:D{2a)n^ 



{2af 



(2.7.26) 



D 



we get the steady state hmit. 



J(oo) = 47rD(2a)noo = STrDaric 



(2.7.27) 



Let us estimate D = 10 '^cw?/sec, and 2a — 10 ®cm, then the time to steady state is 



(2a)! 
D 



— 10 ^sec and is very short. 



Each stationary particle will be hit by, hence coagulate with, SttD a rioo particles per 
second. Since all particles are moving, the steady rate of collision (coagulation) must be 
doubled, i.e., IGnDarioo- As the consequence, the number density of particles, must decrease. 
Each collision reduces the number of particles by 1. Hence 



(in, 



Thus 



or 



drir 



dt 



dt 



-IGiraDn^ where D 



(2.7.28) 



kT 
Qir/ia 



dur 



— IGnaDdt 



which may integrated to 



'^oo(O) 



-IGTraDT 



Note that 



Finally, 



167raD 



IGTra/cT 



Gir/ia 



8kT 
3 II 



^oo(O) 



1 + [16naD]n^{0)t 1 + Konoo{0)t 



(2.7.29) 
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where 



is the coagulation constant while 



is the coagulation time. 

Prom Fuchs, Table 28, p 291. 



8^ 
3 ^ 



1 



coag 



a (cm) 


10"^ 10-6 10-5 10-4 10-3 


Ko X 10^° (cmVsec) 


323 34 5.56 3.19 2.98 



How long does it take for n^{t) to drop to one-tenth of its initial value? 

noo(O) _ 



where 



t = 



K 



KoUooiO) 

_AkT 
3 



It can be estimated that tj_ ~ 125 sec, if 2a = O.lum, and T — 293°ir. 
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2.7.3 Appendix A: Proof of (2.7.19) 

r 1 dn 



dw 

'at 



2a n^o dt 



dw dw/dr ^ (^'^ ^ 
dx dx/dr dr 



1 Ur. 



n 



2a rir. 



r 1 dn 

2a rioo dr 



n r dn 
rino dr 



d'^w 
dx'^ 



2a 



1 dn 1 dn r d'^n 



rioo dr rioo dr n^o dr"^ 



2a 



nr. 



2 dn d'^n 



+ 



r dr dr"^ 



Substituting these results in (2.7.19), we get 

— - D'i2 ? i— -— ^ - — — ir'^ — 
dt \ dr'^ r dr J r'^ dr \ dr ^ 



with 



D = D'{2af 
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2.7.4 Appendix B: Solution of (2.7.19) by the method of similarity 

Let us seek a transformation 

such that the initial-boundary-value problem retains the same form. 

dt dt' dx dx' 

dt ox^ J \ J 

For invariance we require, —2a — —b, a — h/2. Clearly 

^_ X X^/^x x' 



2\/m 2VD'X'''t' 2^fm' 
satifies the requirement. From the boundary conditions, 

x'A" = X^w' = 1 

which requires that 

c = (2.7.34) 

The initial condition as well as the boundary conditon at a;' A" = oo are trivially satisfied. 
The similarity solution is 

( X \ 

w = w{^) = w 



Some algebra: 



2y/m, 



dw I x f 1 \ w'^ 

— w — = w 



hence 



or 



Integrating 



dt dt 2VD' V 2^2/3; 2t 
jj,d'^w ^, w" w" 
'dx^ ~ Im ~ IP 

w% w" 



w" 
w' 

dlogw' _ 
log(w') = -^^ + Const, 



and 



Thus 



so that w{oo) — 0. Since 



The integral is \/7r/2- Hence 



w — ce ^ 



2 dw 



w = —c I e ^ dz 

I z 



J z 

W = 1, ^ = 



=-4 

Jo 



e ^ dz 



c = -1 



w 



1 2 fOO 



e ^ dz 







TT Jo 



